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Abstract. The production of four fermions plus a visible photon in electron-positron collisions is analyzed, 
with particular emphasis on the LEP2 energy range. The study is based on the calculation of exact matrix 
elements, including the effect of fermion masses. In the light of the present measurements performed at 
LEP, triple and quartic anomalous gauge couplings are taken into account. Due to the presence of a visible 
photon in the final state, particular attention is paid to the treatment of higher-order QED corrections. 
Explicit results for integrated cross sections and differential distributions are shown and commented. The 
features of the Monte Carlo program WRAP, used to perform the calculation and available for experimental 
analysis, are described. 
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1 Introduction 

One of the main goals of electroweak physics at LEP2 is 
the study of the properties of the W bosons [|] . The center 
of mass (cm.) energy above the threshold of VF-pair pro- 
duction offers the possibility to extract information about 
the mass and couplings of the W boson from the analysis 
of four-fermion (4/) final states. As well known, because 
of the precision of the experimental measurements, radia- 
tive corrections to e+e~ 4/ processes are needed in or- 
der to provide suitable theoretical predictions ||l|,|| . From 
this point of view, events with four fermions plus a visible 
photon, i.e e^e~ Af + j, are a building block of the 
full 0{a) electroweak corrections to 4/ processes, provid- 
ing the hard bremsstrahlung contribution. Furthermore, 
radiative 4/ processes are also an interesting physics sub- 
ject by themselves, since the luminosity achieved at LEP 
makes them directly accessible to the experimental inves- 
tigation, as recently discussed in refs. where first 
results on the measurement of the W^W^j cross section 
have been reported. 

A very peculiar feature of the processes under con- 
sideration is that they give the opportunity of directly 
testing the non-abelian structure of the gauge boson self- 
interactions. Actually, as other processes studied at LEP, 



such as e"*"e 



W+W~ 



4/, e+e- 



Weve (single-W^ 
production) and e^e~ vv^, e^e~ ^ 4/ -t- 7 reactions 
are sensitive to trilinear gauge couplings (TGC). More im- 
portantly, they can be used to test quartic gauge couplings 
(QGC), since they are, together with vv^^ final states, the 
only accessible LEP2 processes that contain quartic gauge 
boson vertices with at least one photon at tree-level H. 



Quartic vertices involving only massive gauge bosons give 
rise to six-fermion final states and are outside the sensi- 
tivity of LEP, being accessible only at the energies of a 
future e+e" Linear Collider (LC). Both charged current 
(CC) radiative 4/ processes, mediated by two M^-bosons, 
and neutral current (NC) ones, mediated by two Z-bosons, 
are in principle suitable to examine the effect of possibly 
anomalous gauge couplings (AGC). In this work particu- 
lar attention is paid to CC processes, because of the larger 
cross section of WWj events with respect to ZZj final 
states. 

The first tree-level calculations of e+e^ ^ 4/ + 7 pro- 
cesses were performed in refs. In these papers all 
electroweak contributions as well as fermion mass effects 
were accounted for by exploiting different approaches to 
the automatic calculation of the exact tree-level matrix 
element Since then, some event generators for the 
simulation of 4/ -I- 7 events in e^e~ collisions have been 
developed: RacoonWW ||l^, a generator based on the calcu- 
lation of all 4/ + 7 final states in the massless approxima- 
tion ||ll|; CompHEP ^ and grc4f ||l|], which are general- 
purpose packages relying upon the automatic calculation 
of tree-level amplitudes (including fermion masses) and 
phase-space; Helac/Phegas jl^ , a program implementing 
a recursive algorithm for the calculation of the scatter- 
ing amplitudes. The interested reader is referred to ref. 
for more details. In addition to the above computational 
tools, also a calculation of the massive matrix element of 
e~^e~ —> 4/ + 7 processes have recently appeared in the 
literature Pq |, accompanied by a detailed phenomenolog- 
ical analysis of fermion mass effects in 4/ and 4/ + 7 final 
states at LEP2 energies [|l6|. In all the theoretical stud- 
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ies devoted insofar to 4/ + 7 production, the effects of 
quartic anomalous gauge couplings (QAGC), which are 
a window on the mechanism of spontaneous symmetry 
breaking and are presently of special experimental 
interest, have not been taken into account. Actually, re- 
cent phenomenological studies on the subject of QAGC at 
high-energy e~^e~ colliders have been performed by con- 
sidering three-vector boson WWj, ZZ^, Z77 production 
and treating W, Z particles in the on-shcU approxima- 
tion |T8 19 . Anomalous quartic couplings in I'D-fj pro - 



= (27r)6 d<p3{P;p5,Qv,,Qv2)d'^2iQv^;Pl,P2) 
d^2 {Qv2 ; Pz,Pi)dQv^ dQl^ , 



(1) 



duction via WW fusion have been analyzed in ref. 
Experimental searches for QAGC at LEP rely upon the 
theoretical results of refs. ||T^,|l^,pO|j2ll.0 and make use of 



the computational tools of refs. |18lM 

In the light of the present situation and in view of fu- 
ture measurements at LC, a full calculation of e^e~ — > 
4/ + 7 processes, including the effects of AGC and of the 
most important radiative corrections, is desirable. This 
task is accomplished in the present paper, by presenting 
the new event generator WRAP (W Radiative process with 
ALPHA and Pavia) for the simulation of 4/ -I- 7 processes 

at e+e" colliders. The program is based on the calculation = (STr)^ d<p2{P; Qvi , Qv2)d<l>3{Qvi;P5,Pi,P2) 



where Vi and V2 indicate the W gauge bosons, the mo- 
menta Pi with i = 1, . . . , 4 stand for the momenta of the 
final state fermions and P5 is the photon momentum. The 
eleven independent variables have been chosen to be: 

— photon variables E^, 9^ and (jj-y in the cm. frame; 

— invariant mass squared Qy_^ and Qy^ ; 

— three 9 and (f) angle pairs in the rest frame of each 
decaying "particle" , namely in the frames given by the 
conditions P — p5 = 0, Qy^ = and = 0. 

When the photon is emitted from the final state, which 
means for a CC process from a virtual W or from a vir- 
tual charged fcrmion, the following decomposition can be 
conveniently considered: 



of exact matrix elements, including the effect of fermion 
masses, both for CC and NC processes. Charged trilin- 
ear anomalous gauge couplings (TAGC) and the genuinely 
QAGC, i.e. those giving no contributions to trilinear ver- 
tices, are included in the calculation, as well as the large 
effect of initial state radiation (ISR). A tuned compari- 
son between the predictions of a preliminary version of 
WRAP and of the other two event generators RacoonWW and 
Helac/Phegas can be found in ref. 

The paper is organized as follows. In Sect. 2 the main 
features of the the calculation are described. After the 
description of the treatment of the multi-particle phase- 
space, the theoretical details concerning the calculation of 
the exact matrix elements, the implementation of anoma- 
lous gauge couplings and ISR are given. A sample of nu- 
merical results as obtained by means of the Monte Carlo 
WRAP is presented in Sect. 3, paying particular attention to 
the contribution of fermion masses, to the impact of ISR 
and to the effects of AGC at LEP2 and LC energies. Con- 
clusions and possible perspectives are drawn in Sect. 4. 



2 Features of the calculation 

2.1 Phase-space integration 

The kinematics of the 2^5 particles processes has been 
treated generating the 5-body phase-space recursively, since 
the process can be seen as production and subsequent de- 
cay of a pair of massive gauge bosons. 

Concerning CC like processes, the configurations of in- 
terest are related to a photon emitted from the initial state 
(see Fig. |l|), from an intermediate M^-boson (see Fig. 0) 
and from the final state charged fermions (see Fig. 
As far as emission from the initial state is concerned, the 
adopted phase-space decomposition reads as follows: 



d^2 iQv2 ; P3,Pi)dQl^ dQ\ . 



(2) 



In the case of the photon emitted from an internal 
gauge boson the independent variables can be chosen as 
follows: 

— photon variables E^, 9^ and 0^ in the cm. frame; 

— invariant mass squared Qy_^ and Qy^ ; 

— one 9 and (j) W^-angle pair in the cm. frame; 

— two 9 and 4> angle pair for pi and pa in the rest frame 
of the bosons Vi and V2, respectively. 

In the case of photon emission from a final state fermion 
the following independent variables can been adopted: 

— invariant mass squared Qy_^ and Qy^; 

— one 9 and (j) VF-angle pair in the cm. frame; 

— one 9 and (p angle pair for pa in the rest frame of the 
boson V2; 

— energies of pi and p^ momenta in the cm. frame; 

— azimuthal angle (j) of pi in the rest frame of Vi ; 

— in the rest frame of the radiating fermion; 

— cos ^-y-/ in the cm. frame , 

where 9^-f is the relative angle between the radiating 
fermion and the photon. 

An analogous phase-space parameterization has been 
implemented in WRAP for the case of NC processes, neglect- 
ing of course the channels related to the photonic emission 
from internal lines. The above phase-pace decompositions, 
iterated for each possible radiation pattern, give rise to 
several channels, depending on the final state considered. 

In the previous formulas d'Pn represents the element of 
n-body phase-space given by: 



d<?„(P;pi,...,p„) = {2i,)U\P~-Y.p.,)\{ 



d^ 



Pi 



=1 (2^; 



3 2E,- 
(3) 
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Fig. 1. Example of Feynman diagrams for photon radiation from the initial state. 






Fig. 2. Example of Feynman diagrams for photon radiation from W bosons. 





Fig. 3. Example of Feynman diagrams for photon radiation from final-state fermions. 



The code works taking into account all the configu- 
rations discussed above according to a standard multi- 
channel Monte Carlo approach p^ . 

In order to perform an efficient event generation, the 
peaking behaviour of the matrix element has been treated 
in the following way: 

— the squared invariant masses of the massive gauge bosons 
V are sampled according to a Breit-Wigner distribu- 
tion centered around My, while the photon propagator 
is sampled according to the distribution; 

— the infrared divergence is sampled according to a 1 jE-^ 
distribution; 

— the coUinear peak arising from the photon emission 
due to an external charged fermion is taken under con- 
trol by sampling it with a distribution proportional 
to 1/(1 — /3cos6'), where is the separation angle be- 
tween the radiating fermion and the photon, and /3 — 
•\/l — m? / E"^, m and E being the mass and the energy 
of the fermion respectively. 

For a realistic account of gauge boson properties, and 
to avoid integration singularities, it is mandatory to in- 
clude the gauge boson width in the propagators. The so- 
called fixed- width scheme 0,12^ is adopted in WRAP. Actu- 
ally, as shown in ref. ||l^, the fixed- width scheme, even if it 
violates SU{2) gauge invariance, is a reliable U{1) gauge- 
restoring method and is able to guarantee predictions for 
e^e~ ^ 4/ -|- 7 processes in good numerical agreement 
with a scheme preserving all the relevant Ward identities, 
such as the complex-mass scheme |10(|. 



2.2 Tree-level matrix element 

As already mentioned, the present work is based on the 
calculation of the fully massive Born matrix element of 
e"'"e~^4/-|-7 processes. The exact matrix elements for 
CC and NC e+e~ — > 4/ -I- 7 processes are available in 
WRAP. The calculation is performed by using ALPHA Q, an 
iterative algorithm for the automatic evaluation of tree- 
level scattering amplitudes without using Feynman graphs 
(see ref. |^ for a review of the method and of recent 
phenomenological applications). For the processes under 
consideration, a completely numerical approach turns out 
to be particularly convenient not only for the very large 
amount of contributing Feynman diagrams, but also be- 
cause the calculation can be performed in the presence 
of fermion masses without any additional complication. 
This is of special importance for 4/ 7 final states involv- 
ing muons, where the separation angle between muon and 
photon can be realistically set to zero, and a calculation 
taking into account of the finite muon mass is mandatory, 
to avoid coUinear singularities. 



2.3 Anomalous gauge couplings 

Information about the structure of TGC and QGC can be 
obtained by the analysis of 4/ -I- 7 production processes. 
In particular, CC radiative 4/ processes, although char- 
acterized by a lower statistics, are potentially a comple- 
mentary channel to the 4/ final states in order to test the 
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effect of TGC, because of tlie larger amount of diagrams 
involving trilinear gauge interactions. More importantly, 
4/ + 7 processes are interesting in order to put bounds 
on deviations from standard quartic gauge couplings. In 
the following, the theoretical details of the parameteriza- 
tion adopted in order to keep under under control this 
important phenomenological issue are described. 



2.3.1 Trilinear anomalous gauge couplings 

It is possible to take into account the effect of charged 
TGC (anomalous and not) by means of the following la- 
grangian p6| , ^ : 



V 



^Ltgc = 9WWV [ gYV^ ( W-W+'' 

{{^PW-^') - W-" {dPW+'')) V ] 

which represents the most general lagrangian describing 
trilinear WWV gauge interactions, with the exception of 
the operators violating C, P and CP simmetries. Ltgc 
has been implemented in ALPHA, and the presence of anoma- 
lous couplings can be studied, as done at LEP |^^, by 
using the relations [^7|2g|] 



r2 



{Akz - Agf 



= A, (5) 



where Any = Ky — 1 and Ag^ = gf — 1. The Standard 
Model (SM) Lagrangian is recovered for g^ — ky — 1, 
A = 0, gY — 0. Triple anomalous neutral gauge cou- 
plings, considered in ref. |pO| a nd looked for at LEP in 
e+e~ Z"f,ZZ processes are not presently taken 
into account. 



In the above equation the Lorentz structure of the opera- 
tors is given by: 



Wi = a^iFf,,F'"'W+W-P 
W2= a^2F^,.Fi'PW+''W- +h.c. 
Zi= a,iF^,F''''ZpZP 
Zi= a.^F^.F^^PZ'^Zp 
W§ = a^.^F^^Z^'-W+W'P 
= a^,,F^,Z'^PW+''W~ 
Wf = a^.^F^^W+P'^ZPW- 



h.c. 
h.c. 



(7) 



Wi = a^,2F^^W+^^PZ-W- + h.c. 
Wi = a^^sF^^W+PPZpW-'' + h.c. , 

where the are coefficients of dimension M" 2. It is worth 
noticing that, by imposing appropriate relations between 
the flj's, symmetry properties, such as for instance SU (2)c 
custodial symmetry or SU{2) x Ujl) gauge invariance, can 
be guaranteed, as shown in ref. . In the parameteriza- 
tion adopted in ref. the are real coefficients whose 
explicit expression can be directly read off from the corre- 
sponding operator structure of ref. |l9| itself. In particular , 
the coefficients ag and Uc, originally introduced in ref. 
and related to the WW^^ and ZZjj structure, can be 
obtained from the above coefficients by means of the 
following relations 



a-zi 

0-IV2 

az2 = 



" 16cos^ 9^A^ "0 

e „ 
16cos^ e^A'-i "c , 



(8) 



where A represents a scale of new physics. As far as WWZ'j 
vertex is concerned, an additional structure has been pro- 
posed in the literature jlj,^^, whose expression can be 
derived from the above Oi coefficients by means of the 
following relations 



2.3.2 Quartic anomalous gauge couplings 

Quartic gauge couplings involving at least one photon 
are analyzed at LEP g]. In particular, W^W~"fj and 
W~^W~ Zj vertices are probed in WW-f ^ 4/ -I- 7 and 
uujj final states ||,^, while e+e" Z77 processes Q 
are investigated to put bounds on the ZZ^j vertex, which 
is a gauge interaction not predicted by SM at tree level. In 
the present work the operators considered in ref. for 
genuine anomalous quartic couplings containing at least 
one photon, namely W^W~ Z^, W^W~j'j and ZZjj 
vertices, have been implemented in ALPHA, upgrading the 
version used in ref. |g3| for the analysis of QAGC in six- 
fermion final states at the energies of future linear collid- 
ers. The implemented lagrangian include all the relevant 
six-dimensional operators and reads as follows 

Lqgc = W1 + W2 + Zi+ Z2 

wi + wi + wi + wi + wi . (6) 



a-wzc — *i6cose„/i^ ^» 

— _ 

"'"'2=3 — * 16 cos 9^A^ ■ 

On the experimental side, bounds on ao, flc and a„ cou- 
plings are quoted by LEP collaborations via the analysis 
of WW J and 1/1/77 final states [|[^,^ . It is worth notic- 
ing, in passing, that, thanks to the implementation of the 
lagrangian of eq. ^ in the ALPHA code, an improved ver- 
sion of the Monte Carlo generator NUNUGPV is also 
available for the study of QAGC in i^Djj events. 

The result of the implementation in ALPHA has been 
carefully cross-checked by an independent analytical cal- 
culation of all V1V2 V3V4 amplitudes, with Vi = 7, W, Z. 
The check has been performed for all the processes ob- 
tained from W^W~ — > Z7 scattering, by permutating 
particles. 
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2.4 Initial-state radiation 

In order to match the precision of LEP measurements, 
the most important radiative corrections have to be con- 
sidered. Among them, it is wed known that undetected 
initial-state radiation (ISR) plays a major role. It can be 
taken into account in the leading log approximation by us- 
ing the QED Structure Function (SF) approach, in terms 
of coUinear or p^-dependent SF p^ , p5[ . Following re- 
cent work done in ref. p5[ | , both prescriptions are available 
in WRAP, for the reason explained below. When ISR is in- 
cluded via coUinear SF, the QED corrected cross section 
can be written as 



(T 



QED 



(s) — J dxidx2 D{xi, s)D{x2, s) daQ^'^^^ {xiX2s), 

(10) 

by convoluting the tree-level cross section with electron 
SF. However, due to the presence of an observed photon 
in the hard-scattering matrix element, the inclusion of ISR 
needs some care. Actually, since the coUinear SF can be 
viewed as the result of an integration over the angular 
variables of the photon radiation, an overlapping between 
the detected kernel photon and pre-emission photons at 
large angle may occurj^. The consequence is that a double 
counting takes place if higher-order QED corrections are 
naively included by using coUinear SF p5| |. 

On the other hand, it is expected that the bulk of 
the correction is well estimated by coUinear SF, since the 
emission of a photon from an on-shell initial state fermion 
is almost coUinear. However, in order to provide a more 
appropriate treatment of photon corrections and give an 
estimate of the double-counting effect, the SF method can 
be improved by means of the use of p^-dependent SF, 
i.e. by generating angular variables for the ISR photons 
according to l/(p • fc), which is the leading behaviour for 
radiation of momentum k emitted by an on-shell fermion 
of momentum p. In such a scheme, the QED corrected 
cross section can be calculated as 



4/-H7 _ 
'QED — 



da; 1 da; 2 / dc^'^ dc''^^ 



^(a;i, c«; s)^(a;2, 42); s)da4/+i-^, (11) 

where Djx, \ s) is the pj^-dependent SF |34|. According 
to eq. (|ll|), an "equivalent " photon is generated and ac- 
cepted as an ISR contribution only if it satisfies a rejection 
algorithm based on the following requirements: 

— the energy of the "equivalent" photon is below the en- 
ergy threshold for the observed photon, for arbitrary 
angles; or 

— the "equivalent" photon is coUinear to a charged parti- 
cle (i.e. under the minimum separation angle required 
in order to be detected), for arbitrary energies. 

Within the angular acceptance of the observed photon, 
the cross section is computed by means of the exact matrix 

^ The same problem is discussed in detail in ref. for the 
process e'^e~ vu + nj. We refer the reader to ref. for 
more details on the strategy here adopted. 



Table 1. Comparison between WRAP and RacoonWW predictions 
for the massless Born cross section of NC processes at ^/s = 
190 GeV. Input parameters and cuts as in ref. |J10|. 



Cross section (fb) 


WRAP 


RacoonWW 




6.76 ± 0.03 


6.78 ± 0.03 




4.248 ± 0.009 


4.259 ± 0.009 




12.65 ±0.03 


12.70 ±0.04 



element for the considered 4/ + 7 final state. Therefore, 
eq. (^) applies to the signature of four fermions plus an 
isolated hard photon, corrected by the effect of undetected 
soft and/or coUinear radiation. Aiming to obtain a correct 
evaluation of the size of the double counting effects, a 
limit of the present treatment of undetected radiation is 
that only ISR is actually considered. This issue could be 
addressed in a more complete way by using, for example, a 
QED parton shower approach as developed in ref. |37|| , in 
order to describe the radiation from all external charged 
legs, thus including the contribution of undetected final- 
state radiation. 



3 Numerical results 

The aim of the present section is to give some details on 
the technical precision of WRAP and discuss the impact of 
the effects due to fermion masses, ISR and AGC on ob- 
servables of experimental interest. 

In order to test the reliability and the theoretical ac- 
curacy of the event generators, a detailed tuned compar- 
ison between the predictions of WRAP and other available 
programs have been carried out in the context of the four- 
fermion working group of the LEP2 Monte Carlo work- 
shop at CERN H. The comparisons, referred to integrated 
cross sections and differential distributions of several CC 
processes, showed perfect technical agreement. The com- 
parison is here extended to NC processes, as shown in 
Tab. 1^, between the predictions of WRAP and RacoonWW 
with input parameters and cuts as in ref. [ p^ . As can be 
seen, also for NC final states the agreement is excellent. A 
further comparison between the predictions of WRAP and 
those of ref. jl^ is reported in Tab. for several cross 
sections of CC and NC processes, in the presence of finite 
fermion masses and in terms of the same input parame- 
ters and cuts as adopted in ref. ||l6|. Perfect agreement is 
registered for all the considered 4/ + 7 final states. 

The phenomelogical analysis makes use of the following 
input parameters: 



Gf = 1.16637-10- 
Mw = 80.35 GeV 
Fz = 2.49471 GeV 

= 0.10565839 GeV 
rur = 1.55 GeV 



GeV"^ Mz = 91.1867 GeV 

sin2 9^ = 1- M^/M: 
Fw = 2.04277 GeV 
rris = 0.15 GeV 



(12) 

The form used for the propagator of the massive gauge 
bosons is, according to the fixed-width scheme, ~ l/(p^ — 
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Table 2. Comparison between WRAP and the predictions of 
ref. for the massive Born cross section of CC and NC pro- 
cesses at ^/s = 190 GeV. Input parameters and cuts as in 
ref. a. 



Cross section (fb) WRAP 



ref. [1 



ude Ue'y 
udscj 



220.1 ±0.5 
217.5 ±0.4 
78.6 ±0.1 
77.6 ±0.2 
213.0 ±0.1 
208.7 ±0.4 
75.2 ±0.1 
590.0 ±0.6 
5.32 ±0.02 
4.15 ±0.02 
3.175 ± 0.006 



220.3 ±0.7 

218.2 ±0.7 

79.0 ±0.3 
77.5 ± 0.2 
213.8 ±0.3 

209.3 ±0.5 

75.1 ±0.2 
593 ±2 
5.32 ±0.03 
4.18 ±0.02 
3.167 ± 0.007 



Table 3. Comparison between massive and massless Born 
cross section for the final state n'^Vfj,cs ± 7 at ^/s — 200 GeV. 
dn/-f is the minimum separation angle between the photon and 
final state charged fermions; other cuts as in eq. (O). The first 
line refers to the massive case, the second one to the massless 
approximation. 



S-t-f (deg) 


Cross Section (fb) 


5° 


74.294 ± 0.029 




75.732 ± 0.022 


1° 


93.764 ± 0.037 




100.446 ± 0.037 



+ iFM). The processes considered are the radiative 
semi-leptonic final states of the kind e+e^ - 
The cuts adopted are: 



uiqq 7. 



I cos 6*^ I < 0.985; > 1 GeV; 
I cos 6*/ 1 < 0.985; £:/ > 5 GeV; 

> 5°' 

> lOGeV, 



(13) 



where 9-^{i) is the photon(lepton) scattering angle, E^{i) is 
the photon(lepton) energy, O^-f is the angular separation 
between photon and final charged fermions, and Af^j is 

the qq invariant mass. 

Table ^ shows the efi^ect of fermion masses on inte- 
grated cross sections at ^/s = 200 GeV for two different 
photon-fermion separation angles. In the first row, the an- 
gular separation 9^_f between photon and all charged fi- 
nal state fermions is fixed at 5°, while in the second row 
= 1°. As expected, the relative difference between 
the massless and massive cross section increases, going 
from 2% of the first row to the 7% of the second row, 
because of the importance of fermion mass contributions 
when the photon approaches the coUinear region around 
an on-shell charged particle. In the case of a final state 
containing a muon, the separation angle 'd-y-f can be re- 
alistically set to zero, because of different behaviour of 
photons and muons in the experimental apparatus. Ta- 



Table 4. Comparison between massive and massless Born 
cross sections for the final state fi^Vf^cs ± 7 at -^s — 200 GeV. 
d-/-f, with f — q, jj, is the minimum separation angle between 
the photon and final state charged fermions; other cuts as in 
eq. (tsn. In the third column, the first result refers to the mas- 
sive case, and the second one to the massless case. The relative 
difference is shown in the last column. See also ref. [01. 



i?7-9 (deg) 


i?7-M (deg) 


Cross Section (fb) 


H%) 


5° 


1.0° 


90.157 ±0.036 
91.903 ±0.035 


1.92 ±0.08 


5° 


0.1° 


104.777 ± 0.046 
115.004 ±0.044 


9.31 ±0.09 


5° 


0.0° 


105.438 ± 0.045 





ble |j shows the difference between massive and massless 
cross section, with the minimal separation between quarks 
and photon fixed at 5° and progressively relaxing the sep- 
aration cut between muon and photon. It can be seen that 
the massless calculation is still reliable for 1° of minimum 
separation, being the relative difference around 2%, but 
it becomes inadequate when the separation falls at some 
fraction of degree, the relative difference being of the order 
of 10%. Therefore, in particularly stringent experimental 
conditions, only a massive 4/ + 7 calculation can provide 
reliable predictions in the presence of muons in the final 
state. 

Figure ^ shows the line-shape of the cross sections 
of the radiative semi-leptonic processes ^^Vf^ud + 7 and 
e'^i'eud + 7, as a function of the cm. energy in the LEP2 
range. The QED corrected cross section via collinear SF 
for the /x+z/^ud + 7 final state is also plotted. The com- 
parison shows that the contribution due to the additional 



i-channel diagrams present in the e v^ud + 7 final state 
is not particularly relevant for the adopted selection cri- 
teria, small being the differences between the cross sec- 
tions of the two processes. Concerning ISR in the strictly 
collinear approximation, its impact on the cross section is 
at the level of 10 — 15%, which is a phenomenologically 
relevant effect in the light of the LEP experimental ac- 
curacy. It is worth noticing that this result, obtained by 
means of a standard treatment of ISR as tipically adopted 
in the experimental analysis of QAGC in radiative events 
at LEP |3[^, just provides the bulk of the effect due to 
ISR but it is affected, as previously discussed and quanti- 
fied below, by a double counting because of the presence 
of a radiative process as hard-scattering reaction. 

The contribution of initial-state photon radiation is 
also shown in Figs. ||, as a function of the threshold en- 
ergy i*;™™ of the observed photon. It can be noticed that 
the reduction factor due to collinear ISR is around 12- 
13%, almost independent of the photon detection thresh- 
old. However, as previously discussed, ISR in the colHnear 
SF scheme introduces a double counting effect when the 
pre-emission "equivalent" photon enters the phase-space 
region of the kernel photon. In order to get an estimate 
of this overlapping contribution, the comparison of the 
corrections due to the collinear SF and p^-dependent SF 
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is shown. Is is observed that the two prescriptions can 
differ at 5% level for E;™" close to 1-2 GeV, while the be- 
comes smaller and smaller as i?™™ increases. Numerical 
investigation points out that, as expected, the discrepancy 
between collinear and pj^-dependent SF is larger near the 
soft and collinear region and at the level of some per cent, 
thus yielding an estimate of the size of the double-counting 
effect at the level of ISR. Therefore, in the presence of 
particularly stringent experimental constraints sensitive 
to the soft and collinear emission, precise predictions de- 
mand a treatment of ISR able to keep under control the 
transverse degrees of freedom of photon radiation. 

Let us come to the discussion of the effects due to (a 
sample of) AGC. Both integrated cross sections (Figs, g- 

and differential distributions (Figs, ^^[jl^ ) are consid- 
ered. In Figs. M\ the (relative) effect of the TAGC A on 
the e+e^ ^v^ud'^ cross section is examined, by plot- 
ting the relative difference between the cross section in 
the presence of a non-vanishing A coupling and the SM 
cross section (A = ), as a function of the A value at 

= 192 GeV. Fig. || shows a comparison of the effect of 
the A coupling on the radiative fx'^i'^ud+j process and the 
corresponding 4/ final state, which, as already remarked, 
differ for their content of trilinear gauge interactions. The 
numerical results for the 4/ process have been obtained 
by means of the program WWGENPV |Q . For the considered 
A values, the relative contribution is almost the same on 
the two processes (obviously the cross sections are quite 
different), giving a difference at 2 — 3% level only for ex- 
treme A values. Therefore, in the presence of standard cuts 
on the observed photon, trilinear gauge interactions due 
to W radiation in radiative 4/ processes doesn't enhance 
the sensitivity to TAGC with respect to a pure 4/ final 
state. This conclusion is further corroborated by the re- 
sults shown in Fig. where the effect of the A coupling 
is studied for different photon cuts, with the aim of sup- 
pressing the mostly collinear fermion radiation by impos- 
ing more and more severe cuts on the detected photon. 
By comparing the relative deviations shown in Fig. ^ and 
Fig. 0, one can conclude that in radiative 4/ processes W 
radiation can be hardly disentangled from the radiation 
off fermions, being the observed deviations almost at the 
same level for all the set of cuts considered. 

The effect of the QAGC , as defined in ref. |l9|, is 
shown in Figs. for the process fi'^Vfj.udj, as a function 
of the parameter fcg^ at ^/s — 200 GeV. For the scale of 
new physics A, the value A = Mw is used, as convention- 
ally done in the literature. Absolute and relative effects are 
shown in Fig. ^ and Fig. ^, respectively. In terms of the 
coefficients of eq. (0) the coupling can be expressed 



as: 



O-wl 



-e g 



uw 
■ftp 



1^ 



g' cos e„ ]„W 



(14) 



The solid line refers to the complete 4/ -f 7 calculation of 
WRAP with input parameters and photon cuts as used in 
ref. |ll9| and the additional cuts on fermions as given by 
eq. (|13|). In order to compare with the results of ref. |19| , 
the dash-dotted line has been obtained by means of a cal- 



culation of the process e'^e" W^W'j (performed in- 
dependently and in agreement with the results of ref. JT9[ ) , 
by taking into account the suitable branching ratios of the 
W bosons decaying into /x+f^ and ud pairs. The dotted 
line is the prediction as obtained by WRAP, with additional 
cuts on the invariant masses of the two fermionic pairs con- 
strained within 75 GeV and 85 GeV, in order to enhance, 
as much as possible, the contribution of diagrams with two 
resonant W bosons. It can be clearly noticed that, even 
in the presence of cuts on the invariant masses of the de- 
cay products, the complete 4/ -t- 7 calculation differ from 
the prediction of the WW^ approximation, thus proving 
the importance of a full calculation for the extraction of 
meaningful limits on QAGC. 

In Figs. |lo[jl^ the most important photonic distribu- 
tions are displayed using the code as event generator with 
the cuts of eq. ( p^ at a typical LEP2 energy ^/s = 192 GeV. 
In each plot, the SM Born and the QED corrected predic- 
tions are compared with those obtained in the presence of 
AGC. The values used for the anomalous couplings are: 
A = -0.25 and k^ / A^ = 0.01. For the sake of compari- 
son, all the data sample are normalized to the same lumi- 
nosity. The cos 0^ distribution and the distribution of the 
cosine of the angle between the photon and the nearest 
charged particle are shown in Fig. 10 and Fig. |ll|, respec- 
tively. Typical peaking behaviour in the close-to-coUinear 
regions is clearly registered. In such regions, a particularly 
significative impact of the QACG fcg^ is also observed. 
Fig. |l^ and Fig. |l^ refer to the and transverse pho- 
ton momentum p± distribution, respectively, showing the 
characteristic infrared peak. As already noticed in ref. , 
these observables turn out to be particularly sensitive to 
the presence of a QAGC in the region of high energy and 
p_L, being the operator involved of derivative type with 
respect to the photon field. In all the considered distribu- 
tions, ISR introduce sizeable effects if compared with the 
deviations due to anomalous couplings. 

As far as the parameterization of QAGC in terms of 
ao,ac,a„ parameters is concerned, numerical results are 
shown in Tab. | and Tab. ||, at = 200 GeV and = 
500 GeV, respectively. The cross sections in the presence 
of non-vanishing anomalous couplings are compared with 
the pure SM predictions. By looking at Tab. [H] and Tab. |^, 
it can be noticed that the sensitivity of the 4/-I-7 processes 
to QAGC is much higher at the energies of a future LC 
than at LEP2, as a priori expected and already noticed in 
the literature for the WW^ process p^,pj|. 



4 Conclusions 

The production of four fermions plus an additional de- 
tected photon in e+e^ collisions is studied at LEP to test 
electroweak gauge boson couplings and in particular to de- 
rive bounds on QAGC. In order to provide predictions of 
phenomenological interest, an exact calculation of 4/ -I- 7 
processes, including the effect of fermion masses, AGC and 
ISR has been performed. On the basis of the experimen- 
tal accuracy, the contribution of fermion masses and ISR 
has been analyzed in comparison with typical deviations 
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Fig. 11. The cosdjf distribution for the process e'^e ^'^v^ud') at — 192 GeV, where 9.yf is the angle between the 
photon and the nearest charged particle. The samples are the same of Fig. |l^. 
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Table 5. Effect of the QAGC ao/yl^ ac/yl^ a„/yl^ with A = 



Mw, on the cross section of the process e+e 
y/s = 200 GeV. 



tid/i //^7 at 



QAGC 



Cross section (fb) 



Standard Model 
ao/A^ = -0.01 
ao/A^ = +0.01 
ac/A^ = -0.01 
oc/yl^ = +0.01 
fflnM^ = -0.01 

fflnM^ = +0.01 



76.0 ± 0. 
77.0 + 0. 
77.2 ± 0. 
75.5 + 0. 
76.9 + 0. 
76.0 + 0. 
76.0 + 0. 



Table 6. The same as in Tab. || at ^ = 500 GeV. 



QAGC 



Cross section (fb) 



Standard Model 
ao/A^ = -0.001 
ao/A^ = +0.001 
ac/A^ = -0.001 
ac/A^ = +0.001 
a„/A^ = -0.001 
ffln//!^ = +0.001 



25.3 + 0.1 
83.8 + 0.3 
88.0 + 0.2 

41.3 + 0.2 

45.4 + 0.2 
26.4 + 0.1 
26.4 + 0.1 



introduced by AGC. A new Monte-Carlo event generator 
(wrap) has been developed and is available for the simu- 
lation of radiative 4/ events. 

The main conclusions of the present study can be sum- 
marized as follows. The effect of finite fermion masses, as 
analyzed in the fj.'^i'^cs"/ final state, turn out to be very 
sensitive to the separation angle "ff-y-f between photon and 
charged fermions, ranging from about 2% for "d-y-f = 5° 
to about 7% for "d-y-f = 1°. For the realistic situation of a 
vanishing separation angle i^^-^ a massive calculation is 
strictly unavoidable. 

Particular care has been devoted to the inclusion of 
ISR, as a consequence of the presence of an observed pho- 
ton in the final state. The contribution of ISR has been 
studied in terms of coUinear and p± dependent SF. Nu- 
merical results illustrate that ISR introduces corrections 
of the order of 10 — 15% on the integrated cross section. 
However, in order to get a reliable estimate of ISR cor- 
rections and to avoid double counting, p± photon effects 
have to be considered. It has been shown that the double 
counting, affecting the QED corrected cross section via 
collinear SF, may reach the 5% level in a realistic event 
selection and hence it has to be taken into account care- 
fully. A more accurate evaluation of double counting ef- 
fects should however consider also the photonic radiation 
off final state charged fermions. 

Both trilinear and genuinely quartic anomalous gauge 
couplings have been implemented in WRAP, and their effects 
on total cross section as well as on photon distributions 
have been investigated. The impact of TAGC on the con- 
sidered observables does not seem to be very sensitive to 
the cuts imposed on the detected photon, suggesting that 
W radiation is not easily disentangled from the fermion 
radiation. Thus 4/ + 7 final states are not the ideal place 



where to look for TAGC, if compared with 4/ final states, 
which benefit of a higher statistics. On the contrary, these 
radiative processes are significantly affected by QAGC. In 
particular it has been shown that difference in the effect 
of QAGC are present between the predictions of the com- 
plete calculations by means of WRAP and the ones obtained 
in the limit of on-shell W bosons, which is the approxi- 
mation presently used in the literature. A more complete 
investigation of QAGC in radiative events at e~^e~ collid- 
ers is currently in progress. 
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